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1 Introduction 

The main motivation for heavy long-lived particle (LLP) 
searches at the Large Hadron Collider (LHC) arises from 
proposed solutions to the gauge hierarchy problem PQ, 
which typically involve previously unseen particles at the 
TeV mass scale. Hadronising LLPs are anticipated in a 
wide range of physics models that extend the Standard 
Model (SM). For example, these particles appear in both 
A-parity-conserving [2H2UIJ151IS1E11E1IS] and A-parity-violating 
mmm supersymmetry (SUSY) and in universal extra 
dimensions theories mm- 

These particles can be stable, or metastableQthe life¬ 
time of these metastable particles may depend on the mass 
splitting with the lightest SUSY particle or on the size of 
any A-parity-violating coupling IHI- 

LLPs produced at the LHC are expected to be slow 
(/? significantly below 1) and, therefore, should have spe¬ 
cific ionisation higher than any SM particle at high mo¬ 
menta. The ATLAS detector [IB] has a number of sub¬ 
systems able to measure the velocity of charged particles. 
The pixel detector m provides measurements of ionisa¬ 
tion energy loss (dE/dx) whereas the calorimeters and 
the muon spectrometer give a direct measurement of the 
time of flight for particles traversing them. A search for 

1 In this article, particles that do not decay in the active 
detector volume are considered stable, and those that decay 
inside this volume are considered metastable. 


stable LLPs has been performed with the ATLAS detec¬ 
tor [IS] with 4.7 fb _1 of y/s = 7 TeV proton-proton (pp) 
collisions using both the full detector information as well 
as the pixel detector information alone, and has been re¬ 
cently updated with the entire y/s = 8 TeV dataset m , 
but without a pixel-only analysis. The CMS Collabora¬ 
tion has recently published [20] an analysis searching for 
stable LLPs based on the measurement of dE/dx, (3, and 
on muon identification. In CMS, a search for metastable 
LLPs has been carried out by looking for secondary ver¬ 
tices m , or disappearing tracks [223- A displaced vertex 
search performed by the ATLAS Collaboration [23] sets 
limits on metastable particles in a number of scenarios. 
Limits on chargino production, from an analysis search¬ 
ing for disappearing tracks, have also been published by 
the ATLAS Collaboration [Ml . 

The analysis described in this article has sensitivity 
to metastable particles if they have unit charge and their 
track length before decay is more than 45 cm in the ra¬ 
dial direction, so that they can be measured in the first 
few layers of the ATLAS tracker. This measurement does 
not depend on the way the LLP interacts in the dense 
calorimeter material nor, to a first approximation, on the 
LLP decay mode. It can therefore address many differ¬ 
ent models of New Physics, especially those predicting the 
production of metastable heavy particles with O(ns) life¬ 
time at LHC energies, such as mini-split SUSY [25lf26] 
or anomaly-mediated supersymmetry breaking (AMSB) 
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models E3I2H1- A metastable gluino with a mass of 1 TeV 
would be compatible with the measured Higgs boson mass 
according to mini-split SUSY models, which also predict 
squark masses of 10 3 * -10 5 TeV, therefore making the gluino 
the only observable strongly produced SUSY particle at 
LHC energies. In AMSB models, SUSY breaking is caused 
by loop effects and the lightest chargino can be only slightly 
heavier than the lightest neutralino, resulting in a heavy 
charged particle that can be measured before decaying 
into very low energy SM particles and a neutralino. 

Results are presented in the context of SUSY models 
assuming the existence of R-hadrons [29] formed from a 
long-lived coloured sparticle (squark or gluino) and light 
SM quarks or gluons, and in AMSB models for the case 
of long-lived charginos. 

The paper is organised as follows. After a brief de¬ 
scription of the experiment (Section [2]) and of the mea¬ 
surement strategy (Section^, the simulation of the signal 
processes is described (Section [4|. The triggering strategy 
is then summarised and the trigger efficiency is calculated 
for .R-hadrons and charginos (Section |5.1[ ), after which the 
event selection is defined and motivated (Section [5]). The 
data-driven background estimation is then described (Sec¬ 
tion [6]) and the systematic uncertainties are presented and 
discussed (Section [7]). Finally, after a brief description of 
the statistical method used to extract lifetime-dependent 
limits on i?-hadron and chargino production cross sections 
and masses, the results are reported in Section [Sj 

2 ATLAS detector and pixel dE/dx 
measurement 

The ATLAS detecto0 consists of a tracker surrounded 
by a solenoid magnet for measuring the trajectories of 
charged particles, followed by calorimeters for measuring 
the energy of particles that have electromagnetic or strong 
interactions with matter, and a muon spectrometer. The 
muon spectrometer is immersed in a toroidal magnetic 
field and provides tracking for muons, which have typ¬ 
ically passed through the calorimeters. The detector is 
hermetic within its rj acceptance and can therefore mea¬ 
sure the missing transverse momentum (whose magnitude 
is denoted by E™ lss ) associated with each event. A com¬ 
plete description of the ATLAS detector can be found else¬ 
where p]6| . The tracker is made of three detector systems. 
Starting from the solenoid magnet and moving toward the 
beam collision region one finds a «400-thousand-channel 
transition radiation tracker [30] followed by a «6-million- 
channel silicon microstrip detector E2, and finally a «80- 
million-channel pixel detector. The pixel detector is cru¬ 

2 ATLAS uses a right-handed coordinate system with its ori¬ 

gin at the nominal interaction point in the centre of the detec¬ 
tor and the z-axis coinciding with the axis of the beam pipe. 

The ®-axis points from the interaction point to the centre of 
the LHC ring, and the y-axis points upward. Cylindrical co¬ 
ordinates (r, (j >) are used in the transverse plane, /> being the 

azimuthal angle around the beam pipe. The pseudorapidity is 
defined in terms of the polar angle 6 as r/ = — lntan(#/2). 


cial for this measurement and is described in more detail 
below. 

As the innermost sub-detector in ATLAS, the silicon 
pixel detector provides at least three precision measure¬ 
ments for each track in the region |7y| < 2.5 at radial dis¬ 
tances of 5 to 13 cm from the LHC beam line. At nor¬ 
mal incidence, the average charge released by a minimum- 
ionising particle (MIP) in a pixel sensor is ss 20000 e~ and 
the charge threshold is set to 3500±40 e~ for each pixel. 
Signals above this threshold are time-stamped within one 
beam crossing; the hit efficiency under these conditions 
exceeds 99%. When detector data are read out, the time 
over threshold (ToT), i.e. the length of time for which 
the signal is above the threshold, is digitised with 8 bits. 
The ToT is proportional to the ionisation charge [32] and 
its maximum value corresponds to 8.5 times the average 
charge released by a MIP track normal to the silicon de¬ 
tectors and leaving all of its ionisation charge on a single 
pixel. If this value is exceeded, the signal is lost. 

The charge released by a track crossing a layer of the 
pixel detector is rarely contained within just one pixel. 
Neighbouring pixels are thus joined together to form clus¬ 
ters and the charge of a cluster is calculated by summing 
up the charges of all pixels after calibration corrections. 
The specific energy loss ( dE/dx ) is defined as an average 
of the individual cluster ionisation measurements (charge 
collected in the cluster, corrected for the track length in 
the sensor), for the clusters associated with the track. To 
reduce the Landau tails, the average is evaluated after 
having removed the highest dE/dx cluster and amounts 
to 1.24 ± 0.19 MeV/g cm 2 for a MIP [53] . The minimum 
measurable /3y with the dE/dx method is «0.3 for parti¬ 
cles with unit charge and is determined by the ToT over¬ 
flow in any of the pixels in a cluster. 

3 Measurement strategy 

Charged massive LLPs are expected to interact with mat¬ 
ter following the Bethe-Bloch distribution according to 
their /3q. The mass of the LLPs can be obtained by fit¬ 
ting their specific energy loss and momentum to an em¬ 
pirical Bethe-Bloch distribution in the range 0.3 < P'y < 
1.5. This range overlaps with the expected average /3y 
of LLPs produced at the LHC, which decreases from 2.0 
to 0.5 for an increase in particle mass from 100 GeV to 
1600 GeV. The parametric function describing the rela¬ 
tionship between the most probable value of the energy 
loss (dFl/dxMPv) and /3y is [33l : 

dE/dxMPv(P'Y) = ln[l + (\p 2 \hY 5 ] ~ Vt- (1) 

The pi calibration constants have been measured [33] us¬ 
ing low-momentum pions, kaons and protons reconstructed 
in ATLAS and their values are monitored by checking the 
stability of the proton mass measurement as a function of 
time. The calibration is found to be stable at the 1% level 
for all data-taking conditions and detector settings. 

Given a measured value of dU/dx, the mass estimate 
to is obtained from Eq. 0 by numerically solving the 
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Fig. 1 . Ratio of the reconstructed mass, computed as the most 
probable value of a fit to a Landau distribution convolved with 
a Gaussian, to the generated mass, as a function of the gener¬ 
ated mass for stable gluino R-hadrons. The yellow band is the 
half-width at half maximum of the reconstructed mass distri¬ 
bution normalised to the generated mass. 

equation dE/dxMPv(p/m) = dE/dx for the unknown m. 
Figure [T] shows the ratio of the reconstructed mass to the 
generated mass for simulated R-hadrons with masses up to 
1600 GeV. The overestimation of the reconstructed mass 
for heavy particles is due to the pion scattering model 
assumed in the track reconstruction and momentum mea¬ 
surement. A ss3% rescaling of the measured mass is then 
applied in the analysis. The half-width at half maximum 
of the reconstructed mass distribution increases with the 
mass value. This is due to the momentum measurement 
uncertainty dominating the mass resolution above masses 
of 200 GeV. 

The measurement strategy consists of looking for an 
excess of events, compatible with the expected measure¬ 
ment resolution, in the mass distribution for particles that 
are selected as LLP candidates. These should appear as 
high transverse momentum (px) isolated particles with 
large dE/dx. Since there is no trigger based on these ob¬ 
servables, events of interest were selected using the lowest- 
threshold unprescaled calorimetric E™ lss trigger. For sig¬ 
nal events, the E™ lss originates from jets from QCD initial- 
state radiation (ISR) and, whenever relevant, by the LLP 
decays to undetected neutralinos. The trigger efficiency 
for all the LLPs searched for with this analysis is shown 
in Section l5Tl 


4 Simulation of signal 

A number of Monte Carlo (MC) simulated signal samples 
are used in this analysis to determine the expected ef¬ 
ficiencies and to estimate the systematic uncertainties. A 
description of the simulation techniques is presented below 
for both stable and metastable R-hadrons and charginos. 
All simulated events are processed through the Geant4 


[Mj standard ATLAS simulation [25] and digitisation, fol¬ 
lowed by event reconstruction. This includes a realistic 
description of additional pp interactions in the same or 
neighbouring bunch crossings (pile-up). In order to take 
into account residual discrepancies with data, the simu¬ 
lated pile-up distribution is scaled to that observed during 
the 2012 data-taking period. 


4.1 Stable R-hadrons 

Pair production of gluinos with masses between 100 and 
1700 GeV is simulated in Pythia 6.4.27 [22! with the 
AUET2B 137 set of MC tuneable parameters and the 
CTEQ6L1 [22] parton distribution function (PDF) set, in¬ 
corporating dedicated hadronisation routines [22! to pro¬ 
duce final states containing .R-hadrons. Additional sam¬ 
ples of gluinos with some representative mass values are 
generated using MadGraph5 m 1.3.33. Since the Mad- 
Graph samples are generated with an additional outgoing 
parton in the matrix element they provide a more accu¬ 
rate description of ISR and thus a more accurate distri¬ 
bution of the transverse momentum of the gluino-gluino 
system. The gluino samples simulated with PythiaG are 
reweighted to match this gluino-gluino system px dis¬ 
tribution obtained from the MadGraph samples. The 
cross sections are calculated to next-to-leading order in the 
strong coupling constant (NLO), including the resumma¬ 
tion of soft-gluon emission at next-to-leading-logarithmic 
accuracy (NLO+ NLL) [2T1I42II431I441H5] . The nominal cross 
section is calculated assuming a squark mass of 10 TeV. 
The uncertainty is taken from an envelope of cross-section 
predictions using different PDF sets and factorisation and 
renormalisation scales, as described in Ref. [45] , 

Simulated R-hadron events are passed through a full 
detector simulation, where interactions with matter are 
handled by dedicated Geant4 routines based on differ¬ 
ent scattering models. The model described in Refs. [22J 
[47], and hereafter referred to as the generic model, im¬ 
poses few constraints on the allowed stable states. This is 
the only model where doubly charged R-hadrons are pre¬ 
dicted, with a production probability of 0.1%. Hadronic 
scattering is described through a purely phase-space-driven 
approach. A second model, referred to as the Regge model 
in the following, employs a triple-Regge formalism to de¬ 
scribe hadronic scattering and describes R-hadrons con¬ 
taining gluinos according to Ref. [48]. More recent models 
for the hadronic scattering of gluino R-hadrons predict 
that the majority would be electrically neutral after just 
a few hadronic interactions. The third scenario considered 
here belongs to this family and is based on bag-model 
calculations presented in Ref. [25] • This is referred to as 
the intermediate model. The probability for a gluino to 
form a gluon -gluino bound state, based on a colour-octet 
model, is assumed to be 10% [2]. Results are presented for 
the generic (Regge) model for gluino (squark) R-hadrons. 
Variations resulting from the use of a model different from 
the nominal one are taken into account as systematic un¬ 
certainties on the signal efficiency. 
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4.2 Metastable R-hadrons 


The simulation of metastable gluino-based R-hadron sam¬ 
ples is performed in a similar way to that of the stable R- 
hadrons, as described in Section |d.l| The gluinos within 
R-hadrons are required to decay via the radiative process 
g —> gxi or g —> qqxi, using Pythia6. Decays involv¬ 
ing tt pairs are treated separately. Gluino masses between 
400 GeV and 1400 GeV are simulated, with the neutralino 
mass either fixed to 100 GeV or set to m{g) — 100 GeV 
(or m(g) — 480 GeV for tt channels). The gluino lifetime 
is varied from 0.1 ns to 10 ns. 

Results of this search are presented as a function of 
the R-hadron lifetime. Signal MC samples with different 
lifetimes are obtained, starting from those simulated with 
a fixed value of the lifetime, by applying event weights, 
such that the distribution of the proper lifetime in the 
modified sample corresponds to the chosen new value of 
the mean lifetime. The event weight w is given by: 


riRH 

w { t rh ) = IJ 

i 


T 0 

-exp 

trh 




( 2 ) 


where n R H, A) and f,; are the number of R -hadrons in 
the event, the R -hadron mean lifetime as set by the sim¬ 
ulation, and the proper lifetime of the i th R-hadron, re¬ 
spectively. The modified mean lifetime obtained after the 
reweighting procedure is indicated by trh- 


4.3 Stable charginos 

A chargino may be stable in a simplified scenario where 
the mass difference between the chargino and neutralino 
is less than 140 MeV so that the chargino decay to a 
pion and a neutralino is kinematically suppressed. Sam¬ 
ples with long-lived charginos are generated using Her- 
WIG-| —b 2.6.3 [50j along with the UEEE3 |5Tj tune and 
the CTEQ6L1 PDF set. The chargino mass is varied be¬ 
tween 100 GeV and 800 GeV. The chargino is forced to 
remain stable and the other particles in the model are 
set to be too heavy to be produced in s/s = 8 TeV pp 
collisions. Signal cross sections are calculated to NLO us¬ 
ing PROSPIN02 J52]. They are in agreement with the 
NLO+NLL calculations within ~2% |53II54I(55] . The total 
cross section is dominated by direct production of XiXi 
pairs (~ 67%), and by XiXi pahs (~ 30%). The rela¬ 
tive proportion of these two production mechanisms was 
checked and found to be constant at the 1% level over the 
considered chargino mass range. 


4.4 Metastable charginos 

Samples with metastable charginos are produced similarly 
to the samples of stable charginos. The mean lifetime of 
the chargino is set to a given value (r-±), and charginos 

are forced to decay into Xi + 7r± i R the Geant 4 simu¬ 
lation following an exponential decay with lifetime r-± 
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Fig. 2. Efficiency for the calorimetric R™ lss > 80 GeV trigger 
as a function of the R-hadron mass or of the chargino mass. 
Separate curves are also shown for the stable and metastable 
cases. Only statistical uncertainties are shown, which are too 
small to be visible on these graphs. 


in the chargino rest frame. The cliargino-neutralino mass 
splitting is set to 140 MeV. Samples with r-± = 1, 5, 15 
and 30 ns are generated for different mass points. To re¬ 
duce the use of simulation resources, all the samples were 
generated with a jet filter requiring at least one generator- 
level jet with pt > 70 GeV and \r]\ < 5. This choice was 
optimised for the first metastable chargino search ST] . The 
present analysis does not make any explicit requirement 
on the energy of the jets, while it requires missing trans¬ 
verse momentum in the event. There is a strong correlation 
between these two requirements. The residual bias due to 
the jet filter is evaluated and assigned as a systematic un¬ 
certainty. 


5 Candidate selection 

5.1 Trigger 

Events are selected by the E™ lss > 80 GeV trigger, the 
lowest-threshold E™ lss trigger that remained unprescaled 
throughout the 2012 data taking. This is based uniquely 
on the energy deposited in the calorimeters. Figure[2]shows 
the efficiency of this trigger as a function of the R-hadron 
or chargino mass. 

The decay of each LLP to jets and a neutralino oc¬ 
curring within the ATLAS active volume contributes to 
the transverse momentum imbalance and leads to higher 
trigger efficiency in the metastable cases. In the case of 
metastable R-hadrons, the E™ lss also depends on the mass 
of the neutralino produced in the decay itself. The two R- 
hadrons tend to be emitted back-to-back and the same 
holds for the two heavy neutralinos, which therefore ap¬ 
proximately balance E™ lss . If neutralinos are light, high- 
Pt jets coming from R-hadron decay also contribute to 
the transverse momentum inbalance. For this reason, the 
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calorimetric E™ lss is significantly larger when the neu- 
tralino is light. The E™ lss also depends on the lifetime 
of the parent particle as this defines the fraction that de¬ 
cay before the calorimeter and therefore affects E™ lss . If 
the lifetime is shorter than 1 ns, the decay happens very 
close to the primary vertex and the calorimetric E™' ss 
does not depend very much on the lifetime. On the other 
hand, if the lifetime is long enough, the decay may happen 
beyond the calorimeter region and, therefore, the calori¬ 
metric £!™ lss is close to the stable case. The 10 ns lifetime 
is an intermediate case, as decays happen mainly, but not 
exclusively, in the tracker region. 


5.2 Offline selection 

This search is based on a sample of well-measured high- 
px isolated tracks in events with large missing transverse 
momentum. The data sample considered in this analysis 
was collected with tracking detectors, calorimeters, muon 
chambers and magnets fully operational and corresponds 
to a total integrated luminosity of 18.4 fb~ 4 with an un¬ 
certainty of ±2.8% measured using beam separation scans 
following the technique described in Ref. [56j . 

The first step in the selection is the confirmation that 
the event has sufficient E™ lss . The E™ lss variable com¬ 
puted using the offline reconstruction mm, which uses 
refined calorimetric information and includes the contri¬ 
butions of the energy of the muons, must exceed 100 GeV. 
Candidate events are then required to have at least one 
primary vertex with a minimum of five tracks with px > 
0.4 GeV. There must be at least one track associated with 
this primary vertexj^] with at least three pixel hits, mea¬ 
sured over 45 cm in the radial direction, and with trans¬ 
verse momentum px > 80 GeV and r; < 2.5. The set 
of requirements described above, including the trigger re¬ 
quirement, defines the preselection entry in Tables [l] and [2] 

The following additional requirements must be satis¬ 
fied by at least one of the preselected candidate tracks in 
order to select the event. 

The track must be isolated. A track is considered iso¬ 
lated if its distance AR = sj(Aif) 2 + (Ac/)) 2 to any other 
track associated with the primary vertex and with px > 
1 GeV is greater than 0.25. About 70% of these isolated 
tracks are high-px leptons originating from W boson pro¬ 
duction. 

The track must not be identified as an electron [59| 
as LLPs can very rarely (<1%) be identified as electrons. 
The selected tracks are required not to match any recon¬ 
structed electron within AR < 0.01. 

The track must have momentum p > 150 GeV and 
the relative uncertainty on the momentum a v /p < 50%. 
The first requirement improves the signal-to-background 
ratio, while the second ensures good mass resolution. 

3 A track is associated with the primary vertex if its longitu¬ 
dinal and transverse impact parameters are less than 1.5 mm. 
The primary vertex is chosen as that with the highest sum of 
Px associated with it. 


The track must not be a muon originating from a W 
boson decay. Muons cannot be simply identified and re¬ 
jected at this stage, as hypothetical very long-lived par¬ 
ticles would often be mis-identified as muons in the de¬ 
tector. Therefore, to reject muons from a W boson de¬ 
cay, a requirement on transverse mass^] (mx >130 GeV) 
is applied. According to simulation, this requirement re¬ 
duces the fraction of W boson events in the data sample 
to ~ 40%. 

The selected track is required to have specific ionisa¬ 
tion measured by the pixel detector larger than 1.800 — 
0.034|p| ± O.lOlp 2 — 0.029|?7| 3 MeV/gcm 2 . This require¬ 
ment corrects the slight \q\ dependence [60] of the dE/dx 
variable and selects ~ 1.3% of the tracks in the data in¬ 
dependently of the pseudorapidity region. The selection 
cut chosen is the lowest with mass-discriminating power 
(below this, the dE/dx values of all particles are too close 
to the MIP value, irrespective of their masses). 

The above requirements complete the selection for the 
stable particle search. One additional requirement is ap¬ 
plied to improve the sensitivity for the metastable case. 
The highly ionising particles can be matched with recon¬ 
structed jets [61] or muons [62]. Out of 85 candidates, 57 
are geometrically matched to muons (AR < 0.01) and 
26 are AR < 0.07 from a jet. The other two candidates 
have no signals in the calorimeters or muon system in the 
vicinity of the LLP. If the LLPs are stable, they are usu¬ 
ally reconstructed as muons. If the heavy particles are not 
stable, the matching with muons becomes much more rare, 
in particular for particles with a lifetime of O(ns). In the 
search for metastable particles a muon veto is applied, and 
tracks that are matched with a muon are rejected. 

Finally, if more than one track per event passes all 
requirements, the highest-px candidate is chosen, in order 
not to bias the distribution of the variables and to allow 
for proper normalisation in the background estimate. 

Table [l] shows the number of events in data and for an 
example gluino i?-hadron signal for the different selection 
criteria. In data 85 events are selected before and 28 after 
the muon veto. None of the events has more than one 
selected track per event. 

Table [2] shows the yields for the same event selec¬ 
tion as in Table [l] but applied to simulated signal events 
with 1000 GeV gluino R-hadrons that are either stable, 
or otherwise decay to g/qq plus a light neutralino of mass 
m(xi) = 100 GeV, and with a 1 ns lifetime. 

Figure [3] shows the overall signal efficiencies for a rep¬ 
resentative set of simulated signal samples to which the 
full selection procedure is applied. 

When the LLPs decay inside the ATLAS active vol¬ 
ume, Er j? lss increases and trigger and offline E™' ss selection 
becomes more efficient than for the stable case. However, 
as the lifetime decreases, the probability to reconstruct a 
track segment in the silicon detectors decreases dramat¬ 
ically. At a mass of 1000 GeV, these two effects give a 
total efficiency of «15% for the 10 ns lifetime samples and 
«1% for the 1 ns samples, while for stable particles the 
efficiency has intermediate values of ~7%. 

4 m T = x /2p T A“ iss (l - cos(A(/)(E™ iss , track)) 
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Table 1. Observed data event yields at different steps of the 
selection procedure compared with the expected number of 
events for 1000 GeV R-hadrons decaying, with a 10 ns lifetime, 
to g/qq plus a light neutralino of mass m(x l) = 100 GeV. The 
simulated yields are normalised to a total integrated luminos¬ 
ity of 18.4 fb -1 and their statistical uncertainty is also shown. 
See text for details. 


Requirement 

Selected events 
(data) 

Expected events 
(t = 10 ns) 

Preselection 

543692 

112±3 

Isolation 

88431 

102±3 

Electron veto 

60450 

102±3 

High-p 

35684 

91±3 

High-rriT 

6589 

75±2 

Ionisation 

85 

68±2 

Muon veto 

28 

62±2 


Table 2. Expected number of events at different steps of 
the selection procedure for 1000 GeV gluino R-hadrons de¬ 
caying, with a 1 ns lifetime, to g/qq plus a light neutralino of 
mass m(x ?) = 100 GeV, and for stable R-hadrons. The simu¬ 
lated yields are normalised to a total integrated luminosity of 
18.4 fb -1 and their statistical uncertainty is also shown. See 
text for details. 


Requirement 

Expected events 
(t = 1 ns) 

Expected events 
(stable) 

Preselection 

19.5T1.1 

62±2 

Isolation 

9.3T0.7 

55±2 

Electron veto 

9.3±0.7 

54±2 

High-p 

4.4±0.5 

52±2 

High-mT 

3.9±0.4 

43±1 

Ionisation 

3.0T0.4 

37±1 

Muon veto 

2.9T0.4 

- 


>, 

o 
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'o 
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Fig. 3. Total selection efficiencies for some MC samples as a 
function of the R-hadron mass or of the chargino mass. Both 
stable and some metastable cases are shown. Only statistical 
uncertainties are included. 


6 Background estimation 

In order to estimate the background, a data-driven ap¬ 
proach is used. The method uses data to fit the distri¬ 
butions of key variables, taking into account their inter¬ 
dependence, and then to generate a large random sample 
of background events based on the same distributions. The 
choice of the control samples takes into account the mea¬ 
sured correlations between the variables used: p , dE/dx 
and ?y. The dE/dx dependence on the path length in the 
sensor is not linear [60], but depends on 77 , increasing by 
~10% from central to high r/|. The dE/dx also depends 
on the particle /3y via the Bethe-Bloch formula and there¬ 
fore on its momentum, until the Fermi plateau is reached. 
Finally p and ?y depend on event kinematics, as high- 
momentum tracks are more likely to be produced at high 
1771 values. 

Two samples are constructed to describe the distribu¬ 
tions of the key variables. Both selections use the full data 
sample, but with requirements minimising the possible 
contamination by signal events. These control region sam¬ 
ples are the same for the searches for stable and metastable 
particles, except for the rejection of track candidates geo¬ 
metrically matched with spectrometer muons in the latter 
case. 

A first sample ( CR1 ) is selected by applying all the 
selections described in Section [5] except for the high ion¬ 
isation requirement, which is instead inverted to ensure 
orthogonality with the search sample. Otherwise, the kine¬ 
matic properties and overall event characteristics are ex¬ 
pected to be similar to the signal region. With this se¬ 
lection, ~6000 background events are kept in CR1 with 
signal contamination of less than 0 . 6 %. 

A second sample ( CR2 ) is used to obtain the dE/dx 
templates and is selected by inverting the E™ lss require¬ 
ment (i?™ lss <100 GeV), while keeping all the other se¬ 
lection requirements unchanged. This procedure ensures a 
negligible signal content in the sample (signal contamina¬ 
tion of less than 0 . 02 %), and that the selected tracks are 
in the same kinematic ranges of momentum and pseudora¬ 
pidity as the search sample. With this selection, ~ 440000 
background events are kept in CR2 , the majority (~ 96%) 
being matched with muons. The ionisation of all the CR2 
tracks and of the subset with the muon veto applied is 
shown in Figure [4 In their bulk the distributions are 
very similar, thus the larger CR2 sample is also used for 
the metastable search with the muon veto. Nevertheless 
a larger tail, likely due to e + e~ pair production, is seen 
at high dE/dx for the tracks with muon veto, and this 
difference is accounted for in the systematic uncertainty 
evaluation. 

A large background sample consisting of five million 
{ p , 77 , dE/dx } triplets is randomly generated according 
to the following procedure: the momentum is generated 
according to a binned function based on selected tracks 
in CR1 ; the pseudorapidity is generated according to the 
77 -binned (where 77 depends on p) functions based on tracks 
in CR1 and the ionisation is generated according to dE/dx- 
binned (where dE/dx depends on 77 ) functions based on 
all tracks in CR2. 
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Fig. 4. Ionisation distribution of all the CR2 tracks (filled 
circles), and those not matched to a reconstructed muon (open 
squares). The two distributions are normalised to their total 
number of entries. 

For the ~50000 random combinations in which dE/dx 
is larger than the selection requirement, both in the sta¬ 
ble and the metastable scenario, the particle mass m is 
obtained given the {dE/dx,p} generated values, using 
the technique explained in Section [3j The normalisation 
of the generated background distribution to the selected 
data is obtained by scaling the background distribution to 
the data in the low-mass region of the mass distribution 
(40 < m < 160 GeV) where a possible signal has already 
been excluded mm- The normalisation is performed on 
the samples before the ionisation requirement, and its un¬ 
certainty is dominated by the statistical uncertainty in the 
data. The complete procedure, from the key variables de¬ 
scription, to the random generation, to the normalisation, 
is tested on signal-depleted regions. These regions are the 
same as CR1 and CR2 except for requiring tracks with 
100 < p < 150 GeV instead of p > 150 GeV. Applied to 
these validation samples and the stable scenario, the pro¬ 
cedure described above yields a predicted background of 
48.9 ± 0.2 (stat.) ± 0.6 (norm.) events, while the num¬ 
ber of events in the data sample is 49, for m > 160 GeV. 
The same procedure applied to the metastable scenario 
yields a predicted background of 16.9 ± 0.1 (stat.) ± 0.4 
(norm.) events, compared to 20 observed in the data, for 
m > 160 GeV. 


7 Systematic uncertainties 

Systematic uncertainties from several sources affecting the 
background estimate and the signal yield have been eval¬ 
uated. The uncertainties are quoted as the maximum de¬ 
viation from the nominal expectation for the background 
or for the signal in the probed mass range. The actual 
systematic uncertainties are calculated and assigned per 
mass bin. 

The uncertainties in the background estimation can be 
divided into three categories: those related to the partic¬ 


ular choices made for the binning, intervals, and fitting 
functions; those related to the different description of the 
key variables for control samples with different composi¬ 
tions than the search region sample; and those related to 
the stability as a function of pile-up. The uncertainty on 
the background estimate due to each of these sources is 
evaluated by changing the description of the key variables, 
repeating the entire generation procedure, and comparing 
the resulting mass distribution with the nominal one. The 
uncertainties are estimated separately for the searches for 
stable and metastable particles. In each iteration, five mil¬ 
lion events are generated as explained in Sect.[6j The re¬ 
sulting uncertainties are summarised in Table [3j 

Table 3. Summary table of the systematic uncertainties that 
affect the background estimations. All the sources are common 
to searches for stable and metastable particles, unless explic¬ 
itly indicated. The uncertainties depend on the mass, and the 
maximum values are reported. In the limit calculations the ac¬ 
tual value of the uncertainty for a given mass is used. 


Source of uncertainty: 

[%] 

Modification of: 


- Binning in p 

5 

- Binning in g 

5 

- Momentum intervals in g 

2 

- Binning in dE/dx 

4 

- Analytical description of dE/dx 

2 

Different fractions of: 


- Non-leptons in CR2 (metastable only) 

2 

- Tracks with more than 3 pixel clusters in CR2 

1.5 

- Non-leptons in CR1 (stable only) 

4 

- W decays in CR1 and CR2 

3 

Pile-up dependence 

4 

Total background systematic uncertainty: 


Stable particle search 

11 

Metastable particle search 

10 


The uncertainties on the signal yield are summarised 
in Table [2] and can be divided into three categories: those 
on the phenomenological modelling of the signal process 
with Monte Carlo generators; those on the modelling of 
the detector efficiency or calibration; and those affecting 
the overall signal yield. 

The uncertainty on QCD radiation is evaluated for 
stable and metastable f?-hadrons as the difference in ef¬ 
ficiency between the Madgraph and Pythia6 samples. 
This uncertainty is large for the benchmark channels in 
which the E™ lss is dominated by the ISR contribution, 
such as the stable f?-hadrons or the gluino f?-hadrons de¬ 
caying into g/qq plus a heavy neutralino of mass m(\i) = 
m(g) — 100 GeV. The uncertainty on QCD radiation is 
evaluated for metastable charginos using the same proce¬ 
dure as described in Ref. |24| . while for stable charginos 
the same uncertainty used for the stable 17-hadrons is 
used, as the E™ 1SS distributions are very similar. For the R- 
hadrons there is an uncertainty on how they interact with 
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the detector material, and this is evaluated by comparing 
the efficiencies obtained from generating events according 
to the three scattering models described in Section |4T| Ef¬ 
ficiencies for chargino events generated with and without 
the jet filter (see Section 4.4) are compared, and their dif¬ 
ference is accounted for as a systematic uncertainty. A sys¬ 
tematic uncertainty is assigned to the lifetime reweighting 
procedure, and is estimated as the discrepancy between 
the efficiencies obtained for the same reweighted lifetime 
starting from different samples. 


Table 4. Summary table for the sources of systematic un¬ 
certainty considered for 7?-hadrons and charginos. The values 
are separately indicated for metastable and stable cases when 
these are different. The uncertainty depends on the mass and 
on the decay model, and the maximum negative and positive 
values are reported. In the limit calculations the actual value 
of the uncertainty for a given mass is used. 


the uncertainty of the efficiency. Since the pile-up distribu¬ 
tion is different in data and MC simulation, the simulated 
samples are reweighted to match the data, and a system¬ 
atic uncertainty is calculated by varying the weighting fac¬ 
tors. The systematic uncertainty on the pixel ionisation is 
evaluated by comparing the ionisation of simulated and 
real tracks. These tracks are compatible with MIPs, se¬ 
lected with the same requirements as those used for the 
search. Other smaller uncertainties are due to momen¬ 
tum [18] and track efficiency |?P|J parameterisation, and 
electron |65] and muon |66j identification. 

The uncertainty on the signal cross section is calcu¬ 
lated as described in Section [4] The uncertainty ranges 
from 15% (at 100 GeV) to 56% (at 1700 GeV) for 17- 
hadrons, and is ss8.5% for the charginos, slightly depen¬ 
dent on the mass. 


8 Results 


Source of uncertainty 

-[%] 

+ [%] 

i?-hadron 



QCD radiation modelling (stable) 

-28 

28 

QCD radiation modelling (metastable) 

-12 

12 

Scattering models 

-9.9 

6.6 

Lifetime reweighting (metastable) 

-10 

10 

chargino 



QCD radiation modelling (stable) 

-27 

27 

QCD radiation modelling (metastable) 

-21 

30 

Generator jet filter 

-8.3 

0 

Lifetime reweighting (metastable) 

-10 

10 

Trigger efficiency modelling 

-4.5 

4.5 

Ft 1ss scale 

-3.8 

3.5 

Pile-up 

-1.7 

1.7 

Ionisation parameterisation 

-5.8 

0 

Momentum parameterisation (stable) 

-1.0 

1.0 

Momentum parameterisation (metastable) 

-2.0 

2.0 

Track Efficiency parameterisation 

-2.0 

2.0 

Electron identification 

-1.0 

1.0 

Muon identification (metastable only) 

-1.0 

1.0 

Total systematic uncertainty 


Stable R-hadron 

-30 

29 

Metastable R-hadron 

-18 

15 

Stable chargino 

-30 

28 

Metastable chargino 

-24 

31 

Uncertainties on signal yield 



Luminosity 

-2.8 

2.8 

Cross section uncertainty (77-hadron) 

-56 

56 

Cross section uncertainty (chargino) 

-8.5 

8.5 


The systematic uncertainties on the detector modelling 
are dominated by the trigger efficiency modelling, by the 
E!f lss scale, and by the parameterisation of the ionisation. 
Systematic uncertainties related to the trigger are evalu¬ 
ated by varying the threshold and resolution parameters 
of the calorimetric 7?™ lss trigger efficiency modelling curve 
and then looking at the efficiency difference between data 
and MC simulated Z —» g + events. Systematic uncer¬ 
tainties of the 77™ lss measurement are evaluated with the 
methods described in Refs. mm and are propagated to 


The mass distribution is shown in Figure[5]for the selection 
of stable and metastable particles. The data, 85 events 
for the stable selection and 28 events for the metastable 
selection, are compared to the background estimate. In 
addition, mass distributions for examples of gluino and 
chargino signals are shown. 

No evidence of a signal above the expected background 
is observed. The largest deviation has a local p-value of 
4.3% and occurs at a mass of 700 GeV for metastable 
gluino 17-hadrons with lifetime r = 10 ns. 

Upper limits at 95% confidence level (CL) on 7?-hadron 
or chargino production cross sections (and, therefore, given 
a model cross section, a lower limit on their masses) are 
extracted using the CLg method m by scanning the 
signal strength over a suitable range and using the pro¬ 
file likelihood ratio as a test statistic. In the procedure for 
setting the limit, the systematic uncertainties on the sig¬ 
nal and background yields, as evaluated in Section [TJ are 
treated as Gaussian-distributed nuisance parameters. The 
statistical uncertainty on the background distribution also 
takes into account the uncertainty due to the normalisa¬ 
tion. Signal, background and data events are counted in a 
mass window of ±1.4<j around the signal peak, where the 
signal peak and the width are estimated by a Gaussian 
fit to the mass distribution in simulated signal MC sam¬ 
ples. Lower limits on the mass are derived by comparing 
the measured cross-section limits to the lower edge of the 
±1(7 band around the theoretically predicted cross section 
for each process. The resulting lower limits set on the mass 
of the stable and metastable particles are summarised in 
Table [5] for the different searches. 

Figure [6] shows the upper limits on the production 
cross section for the gluino 7?-hadron with a lifetime of 
10 ns decaying into g/qq plus a light neutralino of mass 
m(xi ) = 100 GeV, or a heavy neutralino of mass rn(xi) = 
m{g) — 100 GeV. Figure [7] shows the excluded range of 
lifetimes for the same 77-hadron decays versus the par¬ 
ticle mass. Figure [8] shows the upper limits on the pro¬ 
duction cross section for the gluino 77-hadron with a life¬ 
time of 10 ns decaying into tt plus a light neutralino of 
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Mass [GeV] 



Mass [GeV] 


Table 5. The 95% CL lower limit on the relevant LLP mass 
for the different models considered. Other relevant parameters 
(decay mode, neutralino mass, lifetime if metastable) are also 
shown. 


Particle 

Decay 

m(x l) [GeV] 

t [ns] 

m > [GeV] 

g fi-hadron 

stable 


- 

- 

1115 

b fi-hadron 

stable 


- 

- 

751 

t fi-hadron 

stable 


- 

- 

766 

chargino 

stable 




534 

g fi-hadron 

g/gg 


100 

10 

1185 

g fi-hadron 

g/gg 

m(g) 

- 100 

10 

1099 

g fi-hadron 

tt 


100 

10 

1182 

g fi-hadron 

tt 

™(g) 

-480 

10 

1157 

g fi-hadron 

g/gg 


100 

1.0 

869 

g fi-hadron 

g/gg 

m{g) 

- 100 

1.0 

821 

g fi-hadron 

tt 


100 

1.0 

836 

g fi-hadron 

tt 

™{g) 

-480 

1.0 

836 

chargino 

xi + k ± 

™(xT) 

- 0.14 

1.0 

239 

chargino 

Xi +k ± 

m(xt) 

- 0.14 

15 

482 


Fig. 5. Distribution of the mass of selected candidates, derived 
from the specific ionisation loss, for data, background, and ex¬ 
amples of gluino fi-hadron and chargino signals, for searches for 
stable (top) and metastable (bottom) particles. The expected 
background is shown with its total uncertainty (sum in quadra¬ 
ture of statistical, normalisation and systematic errors). The 
signal distributions are stacked on the expected background, 
and a narrower binning is used for them to allow the signal 
shape to be seen more clearly. The number of signal events is 
that expected according to the theoretical cross sections. For 
both distributions, the bin-per-bin ratio of data to expected 
background is also shown. 


mass m(xi) = 100 GeV, or a heavy neutralino of mass 
TO (Xi) = m (g) ~ 480 GeV. Figure [9] shows the excluded 
range of lifetimes for the same fi-hadron decays as a func¬ 
tion of the particle mass. As shown in Figures [7] and [9] 
the sensitivity of the measurement for fi-hadrons is maxi¬ 
mal for lifetimes around 10 ns. Figure [T0| shows the upper 
limits on the production cross section for charg inos with 
lifetime of 1 ns decaying into x? + 7r ± and Figure [ll] shows 
the excluded range of lifetimes as a function of the particle 
mass for the same chargino decay. 
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m(g) [GeV] 



Fig. 6. Upper limits on the production cross section as a 
function of mass for metastable gluino JGhadrons, with life¬ 
time r = 10 ns, decaying into g/qq plus a light neutralino of 
mass m{x l) = 100 GeV (top) or a heavy neutralino of mass 
m[x ?) = mig) — 100 GeV (bottom). Theoretical values for the 
cross section are shown with their uncertainty. The expected 
upper limit in the background-only case is shown as a solid 
black line, with its ±ler and ±2<r bands, green and yellow, re¬ 
spectively. The observed 95% CL upper limit is shown as a 
solid red line. 




m(g) [GeV] 


Fig. 7. The excluded range of lifetimes as a function of gluino 
mass for gluino R- hadrons decaying into g/qq plus a light neu¬ 
tralino of mass mix ?) = 100 GeV (top) or a heavy neutralino 
of mass mix i) = mig) — 100 GeV (bottom). The expected ex¬ 
clusion, with its experimental ±1 <t band, is given with respect 
to the nominal theoretical cross section. 
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m(g) [GeV] 



Fig. 8. Upper limits on the production cross section as a 
function of mass for metastable gluino J?-hadrons, with life¬ 
time t = 10 ns, decaying into tt plus a light neutralino of 
mass mix l) = 100 GeV (top) or a heavy neutralino of mass 
m[x ?) = m(g) —480 GeV (bottom). Theoretical values for the 
cross section are shown with their uncertainty. The expected 
upper limit in the background-only case is shown as a solid 
black line, with its ±ler and ±2<r bands, green and yellow, re¬ 
spectively. The observed 95% CL upper limit is shown as a 
solid red line. 




Fig. 9. The excluded range of lifetimes as a function of gluino 
mass for gluino i?-hadrons decaying into tt plus a light neu¬ 
tralino of mass mix ?) = 100 GeV (top) or a heavy neutralino 
of mass mix i) = m(p) — 480 GeV (bottom). The expected ex¬ 
clusion, with its experimental ±1 <t band, is given with respect 
to the nominal theoretical cross section. 
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Fig. 10. Upper limits on the production cross section as a 
function of mass for metastable charginos, with lifetime r = 
1.0 ns, decaying into x? +7T ■ Theoretical values for the cross 
section are shown with their uncertainty. The expected upper 
limit in the background-only case is shown as a solid black line, 
with its ±ler and ±2cr bands, green and yellow, respectively. 
The observed 95% CL upper limit is shown as a solid red line. 



m(x;) [GeV] 


Fig. 11. The excluded range of lifetimes as a function of 
chargino mass for charginos decaying to Xi+'k ± . The expected 
exclusion, with its experimental ±1 <t band, is given with re¬ 
spect to the nominal theoretical cross section. 
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9 Summary 2 . 

A search has been performed for stable and metastable 3. 

non-relativistic long-lived charged particles identified through 
their anomalous specific ionisation energy loss in the AT- 4. 
LAS pixel detector. The search uses 18.4 fb ” 1 of pp colli¬ 
sion data at y/s = 8 TeV collected by the ATLAS detector 5. 
at the LHC. In the scenario considered, stable charginos 
with masses smaller than 534 GeV are excluded at 95% 

confidence level, and so are stable i?-hadrons with masses 

’ 7 

smaller than 1115 GeV for gluinos, 751 GeV for bottom 

squarks and 766 GeV for top squarks. For metastable par¬ 
ticles the maximum sensitivity is reached at 10 ns life- g 
time for //-hadrons, where masses below 1185 GeV are 
excluded, and at 15 ns lifetime for charginos, where masses 
up to 482 GeV are excluded. Compared to previous searches 
this search provides the best sensitivity for gluinos with -q 
lifetimes between 3 and 20ns. 
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